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ABSTRACT 


Pa 


The relation of the total structure weight of a 
two-dimensional ring frame to the stress distribution in each of 
the elements is investigated, The total weight cof the two- 
dimensional structures investigated were a minimum when design 
stresses had been attained in all elerents of the structure, 


An expression is developed that approximates the change 
jmmesctress distribution in a merber when additions or deletions cof 
material are made in the member itself and its associated members, 


This expression is used as the basis for a rational 
deésten technique that converges to a ring frame of mintimim total 
weight, 
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Tee LNERODUCTION 


The transverse structure of a wall sided ship consists 
primarily of bottom and deck beams, side frames and stanchions. 

This structure is called the transverse ring frame, The loads this 
structure is subjected to are side and bottom water pressure and the 
ship weights as well as dynamic loadings, 

The structural analysis of this type of structure, in two 
and three dimensions, is discussed in considerable detail in the 
literature, Computer aaron for both the two and three dimensional 
aspects are available. These rapid analysis tools are necessary pre- 
requisites for a rational design synthesis, 

In its most general form the optimization of the transverse 
structure entails the specification of the transverse frame spacing, 
the effect of longitudinal frame spacing, and the distribution of 
material throughout the frame. This thesis deals with the last part 
of the process, the optimium distribution of material throughout the 
ring frame of a wall sided ship once the transverse and longitudinal 
frame spacings are fixed, The optirum structure is defined as the 
least weight corbination of deck beam, side frames, etc., that will 


carry the design loads, 





TI, MINIMUM WEIGHT STRUCTURES 


The end and field roments that a loaded beam element develops 
are directly related to its end restraint, These moment distributions 
are subject to wide variations, The bending moment diagrams for a 
beam element supporting a distributed load are reproduced in Figure 1, 

Full rotational restraint at both ends produces a symmetri- 
cal moment distribution and one in which the maximum moments occur 
at the ends, The field moment dominates when the moment restraint 
Sterne /eaac isezero, The intermediate restraint case is character~ 
ized by a symmetrical moment distribution and identical end and 
field moments. In every case, limiting or design stresses are 
reached only at discrete points in the element, 

The size or weight of any element devends on the maxinum 
moment developed when the design load is eee The wide variation 
in maximum moments is accompanied by corresponding variations in the 
weight per unit length of members. Table 1 lists, for the ideal 
section developed in Appendix C, the variation in weight for the 
structural configurations of Figure 1. The absolute minimum weight 
per unit length occurs when the maximum stress occurs at three points 
simultaneously in the beam element, The specific weight increases 
by thirty percent going from the fixed-fixed beam to the free-free 
or fixed-free beam, Doubling the maximum moment for this ideal 


section increases the section weight by sixty percent, 








The weight ratios of Table 1 underline the fact that the 
total weight of any structural system is Coe to how the individ- 
ual elements interact, A minimum total weight of structure can be 
obtained by substituting members whose moduli satisfy the inter- 
mediate restraint condition, In most cases, the actual moment 
distributions that obtain when the structure is assembled will not 
be identical to those assumed conditions, The final total structural 
weignt will denend on how the material distribution has been altered 
Eowobratmeconvem, vec, “Bhrs, in turn, depends on the design criteria. 
The examples that follow illustrate how the total weight varies when 

| ’ 
the design criteria is that design stress must exist in at least one 
element; the stress in all other members must be equal to or less than 
the design stress, 

Tne first example, Figure 2, contains only elements sub- 
jeecec™to bending. The length to depth ratio cf each element is 
assumed to be large enough to exclude shear effects. 

An initial choice of members based upon the intermediate 
restraint condition and whether the elements carry an external load 
gives: 


2 
1/16 u, L, 


Teen. tie 9 2 
Design 


The structure made up of these members does not possess adequate 
strength, The maximum stress in member two is greater than the 
design stress. The minimum weight members do not form an acceptable 


SEructure, 








This structure can only be brought up to an acceptable 
level of strength by adding material, The weight of the final 
structure satisfying the design criteria will vary considerably de- 
pending on where the material is added, There are two directions 


available, One is to add material so that the moments et the nodes 


ct 


are very small, is will cause the field momant in member two to 
predorinate, The other members will be very small, The second 
approaci is to increase ae size of all the members, This will 
bring about redundant moments at the nodes which are substantial, 
The weight curves of Figure 3 illustrate the trade off 
possibilities that exist when all the members are varied, Large 
increases in the specific weight of members one and three are 
accompanied by relatively small increases in mamber one as the dis- 
tribution factor approacnes one, These effects are summarized in 
tiemtctalsveighivcurve, These curves, for all the eases dooked 2c, 


have definite winizums, When L, = Ly = Sh3, the curve shown in 
Figure 3, the minimum total weight occurred when all member sizes 
were equal and the node momentswere maximum, At this minirum point, 
design stresses were obtained in all members, 

The following example is identical to the preceding except 
that all members have to support an external load, The loading is 


such that the fixed end moment of the middle span is twice that of 


the adjacent members; load symmetry is maintained. 


. 








As before, the structure made up of least weight members, 
Sempespondimg sto the intermediate restraint condition, is not 
satisfactory. All members are overstressed. Sufficient strength 
can only be obtained by adding material, 

Total weight curves for this structure are shown in 
Figure 4, Each of these curves have minimum points which occur at 
k= .5. At this point, all members are of equal size and the node 
moments govern, It is obvious then that the minimum total weight 
occurs when design stress levels exist in all the members, At all 
other values for the distribution factor, one or two members have 
too much material, 


Omly perdine stressesmyere considered 
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The substitution of a vertical member at each of the sunports re- 


On 


quires that direct stresses be accounted for (see Appendix B). Tae 
total weigzt curve for this structure also has a unique minirua, 
This rinicem occurs when design stresses exist in all members, 

The iteration technique developed in this thesis conver- 
ges on a structure whicn, wnen the design load is applied, develops 
design stresses in each member simultaneously, These stresses may 
develop at the ends of the rembers or at any point along the 
members. The structures which this design process converges to 
represents the optinum distribution of material for the given loads 
and the fixed frame spacing, The total weight of these structures 
is a minizum, Any other combination of mewbers will result ina 
higher total weight of structure; the stress in some or all the 


V ~~ Lf o~ _ 1 ne -} oie ~ “~ “4” 
members will be less than the desizva stress. 
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MOMENT DIAGRAMS FOR A BEAM SUPPORTING A DISTRIBUTED LOAD 


FIGURE 1 
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BEAM SYSTEMS FOR EXAMPLES ONE AND TWO 


FIGURE 2 
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IfIL, CONVERGENCE PROCEDURES 


The properties of a two dimensional transverse ring 
frame of a wall sided ship ray be considered similar to those of a 
plane frame. The members of this plane frame are subjected to bend- 
ing and direct stresses; side sway is not considered, In its most 
complex form, this space frame is made up of many bays and decks, 

Picking member sizes that are adequate, i.e. possess 
enough strength to support the applied loads, can be accomplished in 
a straight forward manner, St. Denis [1], presents a rational 
approach to this problem, The refinement of these initial choices, 
i.e. cutting the strength of each member until it just meets the 
design stress, is discussed in this section, 

A two dimensional space frame of more than one deck is a 
statically indeterminent structure, The eee distribution resuit- 
ing from the roe loads depend upon the size of each member of the 
frane,. Changes in these stress distributions are expected if the 
size of one member is altered. The effectsof changing the size of 
a member in a reasonably balanced structure will be greatest in the 
immediate vicinity of the member, It is assumed in the following 
that the stress or moment distribution in a member is affected only 
by changes in members that join the defining nodes of the member. 
This is the same as saying that each member of the frame is the 
same as member 2 of Figure 5. The stress distribution in member 2 


is affected by size changes in any of the members shown, 


a ee 








The refinement problem is this: given a frame with a fixed 
layout and subjected to specified design loads, what distribution 
of material will produce stress distributions such that in each 
memoer the maximum stress equals the desisn stress, This means 
that in each member that has a maximum stress less than design stress, 
GHer section is’too heavy. Material must be taken away or the stress 
distribution altered so that the member carries more of the load, 

If the maximum stress in the member is greater than design stress, 
the member is undersized and the opposite holds, 

An approach to the problem that is suitable to machine 
applications is a selective procedure, After each stress analysis, 
the member in the worst condition, over or under design stress, 
is selected. The size of the member is altered based upon existing 
maximum moment, The structure with the corrected member undergoes 
another stress analyses, This cycle is repeated until all members 
are at the design stress level, 

; A more appealing approach is one in aici changes can be 
made in all the members before each stress analysis, The results 
of each stress analysis are used as the basis for changing unsatis~ 
factory members, 

The implementation of this series approach requires that 
interaction effects be estimated, Unless account is taken of how 
the stress distribution changes when the distribution of material is 


altered, the system will be subjected to undesirable oscillations, 





The member shown in Figure 6 depicts a general member taken 
from the two dimensional plane frame, It carries a distributed 
load, W, and has non-symmetrical end restraint, The stress dis- 
trgeution Bsenolesymmetrical: mode one™ils* assured to be overstressed, 

The information required is how will this moment distribution 
change wien any Or all of the members involved change, At each node 
three conditions have to be considered, an increase, decrease, or no 
enawpe in the joint rigidity. 

An increase in the size of all the members adjacent to node 
two will provide additional rotational restraint. This is accompan- 
ied by an increase in the moment at node two and a decrease in the 
maximum moment, If the changes in the members adjacent to node two 
were sufficiently larse, no change in member one would be required 
due to the altered stress distribution, 

Increasing the size of member one alone will cause the 
maximum moment to increase, If the size of the new member is picked 
based on the existing maximum moment, then it also will be over= 
stressed when the structure is reassembled, 

Table 2 is a summary of the conditions that have to be 
accounted for. Column three has the information needed to make a 
rational choice of new members, An expected increase in the maximum 
moment indicates that a larger than usual correction should be made 
in the already overstressed member. This would be accomplished by 
adding an increment to the existing moment. The new member size would 


Tetlecemeils 2¢ded increment. 


See 





Additional tables would be ere for each of the 
following: a field moment governing, end moment identical and govern- 
ing, end non loaded elements. 

The ilardy Cross mathed of moment distribution offers a 
direct method for obtaining numerical estimates of changes in maxi- 
tun moments due to shifts in material distribution about nodes, The 
expression that is derived below actually embodies the tables dis- 
cussed in the last section. 


A first approximateion to the moment N,, of Figure 5 can 


iz 


be .obtained by carrying out the Hardy Cross method of moment dis- 
@eibutien throughn the first unlocking of all eight joints, It is 
not necessary to require that the joints three throush elght be 
fixed, 

First release the joints adjacent to Satie two and 
distribute the unbalanced moments, U. The carry over from joints 


four, six and eight to node two is; 


Ste, 


CO Mee 2 Ke 42 Ys 


Ue = 27/2 Koo U 


Z 8 


This carry over is added to the unbalanced moment at node two, U.. 


The amount of this summation that carries over to node one after 


node two is unlocked is: 


iy 2 Koy [1/2 (k k 9 U. oe eeu) ee) 


pede | “32 Us * Sapna? Wee! 


The carry over from the joints adjacent to node one is 


Found in the same way. It is; 


EO s=e= 1/2 ks Us - 1/2 k, UL ~- 1/2 k.. U 


at ) 2 amo 
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‘The moment of node one in member 1-2, M 5» can be 


1 


expressed in the following approximate form: 


M = FEA -k West 1/2 k 


12 oy 12 §*5) Us + ky Vs + ky YZ) 


ces: ay, 


ely 2 Koy U, Ly 2 Ei Koy U. -o 174 Koy (kyo U, ae 


Keo U, + Kio U,? 


~ VW4k,, kyo (ke, Ue + kg, Ue +k, U,) + 1/4 ki, k 
1/4 k : lo U 
Js aa | 
(1) 
1G Mio is the maximum moment in the member, the change in 


the moment as the adjoining members and member two itself changes is 


desired. The quantities FEH,. and we depend only on the loads and 


ij 


the geometry and therefore are constants, The expression for M9 


can be looked upon as a function of the distrioution factors, Ki5? 


which in turn depend directly on the section properties, 


ee 1 
M fe (Ki 5 eeoveeveeve k 
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The change in 11 is 
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Applying this to equation (1) results in: 


nN a 1 -- - + t- = Tee ee 
iD Uy, AK, t U/2 (ke US t+ ki, UZ + KS, UL) Ak, 
ae Se ie k 1 oe : 
+ 1/ Ky5 (A 51 U. + A 31 U., a dk) UL) ie U, Ako 
1 Yt 
qe la 2 U, (Ki 5 Ax, + AK), |} S54) ae ales (ke4 U 6 + ky, U, + 
Se Yay) Si 
+ 1/4 Koy (Al Bae Ue, + Ake. Us + Ak,» U,) + 1/4 Uy, 
(Ky 5 Ak a1 + AK 5 Ko) 
(2) 


An adequate 


compa of the orders or magnitude of the 


different parts of this expression requires knowledse of the joint 


unbalance, 


If it is assumed, as it was previously, that the loads 


are reasonably balanced throughout the structure, the structure of 


Figure 5 is again applicable, 


(2) reduces to: 


- Usedk 


At l 


12 a2 


Applied to this structure equation 


_ ng lg 
ih? U., Ako, + 1/2 U,, ( 2 Ak,, + Ak, K54) 
5 } 
+ 1/4 U, (kj, 4k) + Aki, kp) 
(3) 
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This is the same as neglecting the unbalance at joints 
three through eight. The inclusion of these terms in a computer 
progran is straightforward, but of course, would require more 
computational time, This final approximateion accounts for Binces 
in any of the seven members involved. A more convenient form for 


equation (3) is: 


ee = allie = ‘ 1 ic 
ae [-U, @ 1/4 ay) nO oz seca Pa! oie 


- -k 4. 
+ 1/2 [ U, (P- ky) 1/2 Us ky] aks 


(4) 


Equation (4) gives a numerical value for tne expected 


change in My if the adjoining menbers and member itself change size, 


j 
This is an approximation as to how the moment distribution in member 
ij is affected when any or all the members change size, This is a 
valid approximation for both loaded and unloaded elements, There 
is no restriction to members undergoing only bending stress; it is 
equally useful for members with combination stresses, If different 
types of materials are used, their effects are included in the 
distribution factors, 

The use of this expression to find changes in the field 


moment is a little more involved, [Equation (4) has to be evaluated 


for each end of the member in question if the moment distribution 





is not symmetrical and added to the existing end moments, These 
new moments combined with the existing end shears give an approxi- 
mation to the change in field moment, 

A word about signs is appropriate, Clockwise moments 
are positive, counterclockwise moments are taken to be negative. 
The signs that result when ire is calculated cannot be interpreted 
literally, A positive value for sai may indicate an expected de-’ 


@rease in the moment, M It is necessary to check the sien of 


se 


the existing moment and then interpret the sign of eee 


ah ee 








A SUBSECTION OF THE TWO-DIIEENSTONAL RING FRAME 


FIGURE 5 


tee 








eee TT 117 


A GENERAL BEAM? HEMBER OF THE TWO-DIMENSIONAL RING FRAME 


FIGURE 6 
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TABLE 2 


Gupmiative frtects oo. Mfember Variations 


(2) 
CP) Member (3) 
Adjacent to be Exp. Change 
Hemoers Chanced inetlex, Ton, 
- aioli 3 a 
End Increase 
w/ + + Ne Chance 
Mom)... i 
Min 
Oo + Increase 
= + Decrease | 
End 
ws - + No Chance 
Mon) 
_ O + Decrease 


Explanation of Table 
This table is drawn up for Member Two of Figure Five, It is 
assumed that erreuititent distribution is nonsymmetrical and that 
Member Two is overstressed, Column three gives the expected change 
in the maxinum moment based upon the cumulative effects of member 
changes summarized in columns one and two, 
Column 1; The sign indicates how the overall stiffness of 
the joint is expected to change based upon the 
Stress conditions in the adjacent members, 
Column 2: The size of Member Two is expected to increase 
Column 3: The expected change in the maximum moment based 


upon the data of columns one and two, 
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IV. TUE DESIGN SPIRA 


The result of the previous section is a technique for 
estimating moment changes. These roment chances are brought about 


when the members at a node change their size, 


65) 


This is of use in the design process because it presents 
a method of estimating not only the magnitude of changes in the 
stress distribution (maximum moment), but also the direction of 
the changes, It permits a more rapid convergence to the final 
structure because each member of the structure can be modified (if 
necessary) on each pass, And also because the members are picked 
based upon conditions that will prevail; not those of the last 
stress analysis, It is essentially a method of predicting the 
final moment ane eaeen and picking the members based on this 
distribution, If the moment changed predicted were exact instea 
of approximations there wouldn't be any iterative procedure 
necessary, 

A desirable aspect of this technique is that it is self 
modifying when applied repeatedly, The predicted moment changes 
are apdeiedl 26h time a new member is picked. In this way oscilla- 
tions are avoided and smooth convergence is assured, It is not 
necessary to perform a new structural analysis until this self 


modifying aspect is completed, This occurs when there are no 


ne 
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significant’ differences in the predicted moment ae as from one 
¢ycle to the next. | . 

The rational design process developed in the sections that 
follow converges to the optimum distribution of material for a trans- 
verse ring frame, The general procedure followed is to use the 
results of a stress analysis as the basis for modifying the material 
distribution in such a manner that desien stresses will be attained 
in eaci member, The results of Section II indicate that a structure 
which develops design stresses in each member simultaneously is 
the least weignt structure, 

- 

Equation (4) is used in this process to make corrections 
to the results of the stress analysis. <A change in the ecantlines 
of any nemoer of the ring frame changes the stress distributions, 

As pointed out in Section IIt and summarized in Table IL, the 
ettect of a change in scantlings depends not only on the magnitude 
of the change, but also the existing moment distribution, Once the 
maximum moment in the member is identified, equation (4) gives 
direct information as to how this moment will change when tne 
scantlings of the members of the ring frame are revised, 

This correction to the existing moment distributions 
provides the means of facilitating rapid convergence to the optitun 
structure, Chanses in member scantlings are based on these pre- 


dicted moments, Tne objective of design stress in each member of 


the ring frare is obtained much more readily. 








-The use of equation (4) in the design process requires 
the standard output of a stress analysis, end forces, plus end 
moments, The process begins with the results of the stress analysis 
which used the initial choices of members. 

The computer begins aren member one, The governing 
moment is ascertained and the stress level is computed based on this 
maxirum monent, If tne stress level is within a certain prescribed 
range of the design stress, the member is not changed, This pro- 
vides a means for ending the iterations. If the stress level is 
not satisfactory, a new giver size is calculated and stored in the 
matrix, K, This procedure is repeated for each member of “the 
structure. The result of this cycle is the matrix k. The members 
that make up K represent the first cut at bringing about the desired 
distribution of material, These merbers are fictitious in the 
sense that they are not used in an actual stress analysis, The 
properties of these members are based upon tne existing stress 
distributions, 

The computer now starts to revise the members of matrix k;} 
it starts with the first loaded member, The stiffnesses, ae 
the members at the defining nodes of the member to be changed are 
computed using the data of matrix K, This data is used to compute 
ae 


field moment governs Au 


or ep or both depending on the stress distribution, If the 


13 is computed for each node; when an end 


moment governs OG is calculated for the governing node, 





The maximum moment of the member to be changed is modified 


by adding the correction factor AM,,. The properties of a new 


al 


Cu. 


section are calculated based on this modified moment, These 
properties are inserted in the proner slot of matrix k, 

The computer chooses the next loaded member anc repeats 
the above, The unloaded members are revised after new properties 
have been picked for all the loaded members, This requires fewer 
steps because the end moments always govern, 

The result after all the members have been looked at is 

? 
a new. matrix of menbers, Kk , A stress analysis could be initiated 
a 
at this point. The results of the examples indicate, however, that 
a smoother convergence is cbtained if this process is repeated, i.e, 
" 

fond the m@agitrixihie. 

A stress analysis is now undertaken using the nenbers of 

4 , 

the matrix K , The results of this stress analysis are used as 
the starting point for a new design cycle, 

Flow charts containing the details of the design spiral 


are contabmed im Appendix A, Illustrations of the design spiral 


in action are contained in Appendix 5B. 





V¥. RESULTS AND CONCLUSIONS 


A rational design procecure has been developed which 
permits tne@raped syitnes#s of tWemtrausverse ring frame of a 
wall sided ship. The resulting transverse structure reoresents 
the optimum distribution of material, the minimum wieght combus- 
t#0n of deck beams, side frames, etc., within the constraints of 
fixed longitudinal and transverse frame spacing, 

Detailed applications of this procedure are carried out 
in Appendix Bb, Each step follows the subroutines develoved for 

, 
computer adaptation of this technique, No attempt has been made 
to adapt the procedure to hand calculations, 
This method converses rapidly, Large member changes were 
required in each of the examples to attain the design stress in each 
member, Example I required only one additional stress analysis 
after the initial choice of members were analyzed; Example II re- 
quired two additional moment distributions, 

The speed with which a satisfactory ring frame is obtained 
depends directly on the following: the limits set ee which the 
member stress must fall and the use of weighting functions in con- 
junction with the predicted moment changes, 


The limits set for member stresses in Appendix B were 


liberal, This, of course, was conducive to a rapid convergence, 





These limits are not fixed, however, and can be set at any desired 
level] ,~eit is™pessibile to puitetew limits on the members stress ‘with 
each problen, 

The predicted moment changes in the examples were used 
as given by equation (4), The true approximation in eauation (4) is 
in the magnitude of expected moment change, The direction of the 
change, plus or minus, is exact if the change is greater than a few 
percent, In all the cases investigated, the magnitude of the change 
was always low, Convergence would have been more rapid if the 
changes were multiplied by a weighting factor, This indicates that 

- 
weighting factors, if necessary, can be developed as experience is 
gained with the basic program, These weighting factors would become 
part of the basic program, 

The results of Example II of Appendix B indicates that 
members whose continuity is not required to carry the applied loads 
may be eliminated, This is easily avoided by applying nominal 
loads to these members, e.¢, their own weight, In an actual 
structure, this will seldom be a problem, 

The use of the ideal section of Appendix C was a necessity 
in the development stases of this program, The implementation of 
this design procedure, i.e, the actual programming, does not re- 
quire this ideal section, If a tape of actual sections is available, 


the program can use it and still work, The effect tnis would have 


es 





on the sveed of convergence would be mixed, The actual members in 
most cases have larger section moduli than those the program would 
pick, The use of actual sections requires,in addition, the develop- 


memtmemea trial and @@ror soliteion in the subroutine COMSTRESS, 


ee ONG es 





APPENDIX A 


An algorithm covering the complete design process is 
developed, Flow charts covering the details of the individual steps 


of the algorithm are also presented, 


Design Algorithm 


Start with member one compare the maximum 
ana desicn steesces, lf the maxirem stress 4s 
within the prescribed range, insert the pre- 
sent merber properties in matrix K, If maxi- 
hum stress unset... picm™aiewerember, Insert 


its properties in the matrix k, 












Repeat the above for each merber, If all 
members have satisfactory stress levels the con- 
puter stops - the design is completed, If any 
or all members were changed, the iteration 
process is begun by returning to the loaded 


member number one, 


Required stiffness factors are calculated 


using elements in matrix K, The moment change 


15 détemwmemedewaune this datay The existirs, 


maxirum moment plus this morent change is the 


® 
it Cas Sen = ba - : ap ora 2 ox we 7 = 
new em? ee Bic we Lhe os ey et .? “= 


‘ 


this desien mowent and the direct stress. 





——— 









The new member is placed in matrix kK, 

















This process is repeated for all members 
even if stress levels are initially satis- 


factoum  itheseatrimecontaining these new 


' 
members is now designated k , 


ee ee ee ee en ee eee eee 
SE Ta 


i] 
The data of matrix K is now wsed to 


calculate the required stiffness factors. 
The process is repeated, The matrix 


vv} 
formed is Kk 





A stress analysis is carried out using 
Giemeer ~eeime™reenits of this analysis 


are tne starting point for a new desipfn 


Cre ve > 


The flow charts are made up in the form of subroutines, 
They cover all the steps presented in the design algorithm, 

Tree trateelow cart 15 subroutine torent. It uses as an 
input the data of the preceding stress analysis. This subroutine 
primarily finds the maximum moment and calculates tue maximum 
stress, When field moments govern tne member is given a special tag, 

Subroutine Change checks the maximum stress in each member, 


It stops the iteration process if all members are satisfactory, 


- 39 - 





Subroutine Stiffness uses the data of matrix K and the 
data of initial member choices to calculate the expected changes in 
stiffness, It also determines the moment change and adds this to 
the maximam moment, The new member parareters are obtainad by 
calling the apprepriate subroutine, t,e, Pick er Comstress, This 
subroutine assumes that a member numbering system is available which 
facilitates picking the required data from the matrix K, etc, 

Subroutines Pick and Comstress both choose new member 
properties, Comstress handles members that have significant direct 
stress as well as bendine stresses, 

The calling sequence for these subroutines would be: 

CALI, SUBROUTINE MOMENT 

CALL VOUBROUTINE CHANGE 

CALL SUBROUTINE STIABNESS 


CALL STRESS = or any suitable stress analysis prograr, 





SUBROUTINE MOMENT. 





N=—-Total # of Members 


Set Member Counter, n = 1 





4 > a . 
= CE Stop, Transfer 


Control to Main Program 


— 





SUBROUTINE CHANGE 


N<«— Number of Members 


n<—— | 












Puts Present Memb 
Properties in 






er 
be If Member Flagged for 
Combined Stress, Sub- 
POUtine Comeeress 


Called, Otherwise 









SM@broutine Pic: Called, | 








Seen 





DECo Erimt out 
Matrix kK. 
Problem Completed 


Stop. Transfer 


COntscol «o lain 





SS te a Otte uy” 





SUeOUTINE STLFFNESS 


N <— Number of Members 


n<«—1l1; F <—l 















Same as Left Side 
Pecept Dittereac 


Members, 





Repeat the above for 
Node j 





Finds M..)., 
a4 ria 


Subroutine Comstress 
Called if Nember rlaggec. 


Otherwise Subroutine Pick 


Control Ca, shee. | 
Analysis Proaren, 


eg ee 





SUBROUTINE PICK 










3 
Ds, ——— (25/032) °| 


Rin pens OD 2) a 


1J 





ab 


Return Control to 
Calling Program 


SU ROU COMMS TRESS 


ee ee os oe 





a <—- F,,/o x ,052 
ij 6 






SE A 
Y <—- ch i a a 
Ge 032 
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, tae » 032 p? 
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fieous.er Control to 
Calling Subroutzne 
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APPENDIX 


EXAMPLES 


' Gere iar 


EXAMPLE 1 


ae Bowe 


The seven membered frame shown on rage 45 is used to 
illustrate the optimization technique developed in the body of 


Ghie tnesis, 


oe ecmlen Oo: tatttal Member Properties 


Cates 


" The initial set of member properties is part of the data 
input to the computer, The following were selected for this 


exanole, 


(a) , 240 ee : I, = 20090 om 


N 
i 


(Hoe Zoe =ee =. 120 ae eee = 76703 
Cee fee = 2. = 0 in : L, = 820 in 


Pvece preyeerrces vere Opeained Using Ere ideal section of Appendix C. 
Bu Stress Aqnivsis 

The Hardy Cross method of moment distribution was used to 
obtain end moments and reactions, The results of this first enaly- 


sis are: 


he = 
ho Moy 


Minx 7 684 KIP FT 


= »% = \N = r A 
ea em 2 = 495 KIP PT 


ae a 


= '. = é 
1B Cami = 180 KIP FT 


fee 4 





C. Subroutine Change 
(1) The maximum.stress in each member is compared to 
the design stress = 40 KSI 


(2) The size of each unsatisfactory member is changed, 


The maximum stress in each member is: 


= hone 
oO 34° koe 
Ca) ax = 50 KSI 
ee = 15 KSI + 10 KSI = 25 KSI = a combined stress member, 


The new properties for each member are: 


ML), 
Zee = 04 in; 1 = 1600 in 
Z ej : 2 
Design 
zZ, = 148 men I, = 1100 in’ 
Z, = 67 ine = 6,3 Rec tT, = 415 ae 
eae |e aif 


The revraining member sizes follow by symmetry, These preper- 
ties make up the matrix Kk, The subroutine PICK would ba used 
to obtain the properties of members one through three; sub- 
routine COMSTRESS would be called to obtain the properties 
for members four through eight. | 

D,__ Subroutine Stiffness 


t 
(1) Calculates mewestifiness factors, ae 


(2) Calculates the change in stiffness factors with 


respect to the initial member choices, 


eee 





D, Subroutine Stiffness (Cont'd) 


(3) Determines 25s and the new factiticus maximum moment 


(4) Picks new members based one the new maximum moment, 


This subroutine begins with the first loaded member, member 1 


t 
ana Calewlates k 


i 

es a8 10 a 
kj, = 184 Lae 0 
dk), = + «126 dk), = 0 


These results are used in equation (4) to obtain AM)» 


A =+ 50.4 K' : M4. _ = + 495 + 50.4 = 545 KIP ET 


M 
“1-7 1-7 


Subroutine Pick now comes up with 


? ? 
Z, = 153 ce I, = 1140 ae 


These member properties are imrediately inserted in the 
matrix K, This provides a constant updating of the data, 

The computer proceeds thru the structure doing the 
above for each member 


Member 2 


: ' 


Kio = +448 ky, = 0452 
Ky. = +450 kag 2490 
t t 
dk, = ~.002 dk, = + .002 
' t 
AM = 1,1K'; - M,. = 685 KIP FT 


{=-2 


v t 


4 
Z = 204 oe: fj = 2000 in 


- 39 = 











D. Subroutine Stiffness (Cont'd) 
Moiber 3. 
' ' 
Kyi = 310 Koy = 0 
- = 5 / . =: 
Ki Loe Koy 0 
te = = 
AK, 5 126 dk) 0 
i] ? 
At = + 50.4 KIP FT MU = 545 KIP FI 
’ ? , 
Z. abel. ‘ce IL = 1140 ant 
Hember 4 
' ? 
Ki 3 args OL? key = 0 
Ki3 ee 4 Ke = 0 
' ? 
Ak} 3 bea ~067 Ake 4 = 0 
? ? 
27 TES gin > Nee 
? 2 ? 4 
Z., = 255 (ine I, = 310 in Subroutine se@ro RESS 
Nemvenes 
eee se t 
Kis Sibyl Key = 0 
Ak) = ,184 Ake) = 0 
? 
Ak; = -.063 Key = 0 
? ? 
bit. = eee LP ET Me =) ee 
rs . 3 aa ee ee yall 4 
Z. a ein, Ae aye Sins T. = 310 in 
fembers 6 and 7 follow by symmetry. 


aoe gee 





: ' 
pie rer a heat this point becomes the matrix K , 


Subroutine Stiffness repeats the above process and produces matrix 
ib 
K ., The results of this procedure are summarized below 


Menber z Gees a (in?) L (in’} 
iL 167 acd ati 1300 
2 204 ahaha 1600 
4 50 6.9 280 
5 50 Ga? 280 


iT 
The matrix kK is now the matrix of actual member 


properties, A stress analysis of the structure made up of these 


members is performea, The results of this analvsis are; 


i, = Deoenele PT G.U=—42°nol = 107 te 
Mo = GeO CLP FT o, = 40 KSI = o 

2 2 O 
MW, @& .41 KIP FT Gye SOvhol = 00 

4 4 Oo 


Stresses and moments for other members follow by symmetry, 
Suvroutine Change now finds all the stress levels 
ty 
satisfactory. The members of matrix K represent the optirum 
structure, The desipga spiral is complete, 
EXAMPLE 2 
This example was chosen to illustrate the procedure for 


a nonsynmetrical system and one in which the field moment governs, 


at least initially. 


ee Oise 





The frame analyzed is shown on page 46. The optimization 
process proceeds as illustrated in detail in example one, 
Bs eee les 


4 
Z, = 240 ieee 2000 in 


All other members - Z = 40 Figg A = 6 ree P= 210 ee 


oe stress Agalvysis 


rs 
i! 


"160 KIB eG? 


7 oO Oly eee 


t aed 
am) 
f 


M, = 168 KIP FT 
M. = 84 KIP FT 


Subroutine Change 


$ 
Pe ee epee 





q 3 | 4 
= ae = i ° = 
oy 50 KSI zy 5G) mah gues I, 2902 
Oo, = 39° hs, te governing : 3 4 
Zz, = 240 in ; I, = 20 Orin 
oe 50 wSI | 
Combined Stress ; 3 ; , : i 
eos Se) eens zy, com lia, A, eNO gene c: I, = 569 in 
q ; 3 - 9 = : ae q z - ds 
zZ« 3 te A, = 5,6 in; I, = 190s 
D. Subroutine Stiffness 
¢ 
Matra. 
] 3 i] a ; ] 4 
Member tain} No Gla): L wean 
ii Doe L Saas 520 
2 208 eo 1700 
4 140 13,5 1050 
a Zz 4,1 100 


ary ate 





The field moment governed in member 2 at this point, 


However once the second set of members were obtained, the end 
moments predominated, The results of the second pass is 


matrix k 
Member 2 ne) A Gre) a Gia 
1 206 anaes 320 
2 203 ———— 1650 
4 190 16.5 1500 
5 1 3 _ 60 


| oan stress Analysis 


Hee oO heme I 
= 655 KIP FT 
M, = 114 KIP FT 
i = 70 


Fe Subroutine Change 


Te hss 


6. = 3876 7551 = ,97 

1 O 
o. = 30 KSI = .75 

2 O 
G = Go nor = 1.25 

b O 
o,. = Q 


On Hy ee 





These stress values require another pass, 
obvious, I 
The second iteration with subroutine Stiffness 


followings watrices. 


Matrix K - second pass 
Member 2 Gin A Ga) 
i 17 ae 
Zz 203 =e 
4, 262 21 
5 = on a me om oe me oe 
" 
Matrix K - second pass 
Member 2 (in?) ne aS) 
= === -—-= 
Z 192 -—-= 
, 267 Zi 
5 tai ae, 
G. Stress Analysis 
NE Pe 1S TID Tt 
M, = 646 KIP FT 
H, Subroutine Change 
oe Soo sol 6 CUS oh 
aa Pome tl.) = 40,7 KSI = 1,025 a 


mons ae 


Members two, four and six are taking 


iiawee trend is 


a1k the load, 


produces the 


iL a) 
40 
1650 


2000 


= oe os me 


aL Boe) 
1550 


2000 


wa oe & 2 





These stresses are within the specified limits, 


fe 
The design spiral is terminated, The members of inatrix K , 


second pass, form the optimum structure, 








3 a 


ol 5 ede a 
eae FI 6 All L’s equal, 


tas 20) 
12 eS yee TZ Le s/t 
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SERUIGHURE FOR EXATIPLE Bel 


FIGURE B-l 
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SHRUCTURE FOR EXAMPLE B-2 


FIGURE B-2 
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AEE EMNDIX C 


THE IDEAL SECTION 


The section moduli of available steel sections listed in 
the Manual of Steel Construction do not vary in a continuous 
manner, The use of the least weight sections compounds the 
problem, Each section is useful for a wide range of section moduli, 

These facts make it very difficult, if one is Wen Ge 
sections listed in the AISC Manual, to pick out trends, Also the 
fact that there are no mathematical relationships existing between 
the important properties of the sections compounds the problems 
encountered when using actual sections, 

A mathematical section, i.e. one in which there are 
specific relationships existing between the section parameters, 
avoids these problems. The section see Sane for this thesis is 
outlined below, It conforms well to minimum weight sections listed 
in the AISC Manual for larger section moduli. At the lower end of 


the scale there 1S a greater divergence from actual sections, 


Let ae flange thickness 


ae = web thickness 
D = depth of section 
We = flange width 


The relationsnips chosen for the section are; 


te = 2 t. 

DB = 50 ¢t 
W 

We = 20 t 


ee fy eee 





The expressions for the area, section modulus and moment 


OF inertia are: 


A= 62 D- 
Z = ,032 D> 
i ae Ae 


A log-log plot of these expressions is made on page 49, 


ee 





200 


(iN3) 


100 


30 


25 





2000 
I 
(1x*) 
1000 
. ; ae 
500 s 
r/ 
Me 
= 
»/f 5 
i a 
f hey 
250 " 
4 «5 10 


BEAM DEPTH (INCHES) 


FIGURE Cel 


ae One 


a 


ly 


Vina 25S 


20 


200 


(ey a) 
100 


25 


510) 





APPENDIX D 


BIBLIOGRAPHY 


Deo ds ee eee wo emmcuure) Desien Gt tre Midshin Section, 
BID Reser C Ceg55,70CE. 1954. 


Bruhn, J., "The Transverse Strength of Ships", Trans, I. N. A., 
Tolan Gene 05257 (lool), 


Hay, W. L., "Some Notes on Ship Structures’, Trans, I. N. A., 
vol, 87, pp. 81-94 (1945). 


DRUestonge Js Ns Gy, ransverse Strength, Trans, T.N. A., 
Vol~enuOre po, 40/=41/ (Oct, 1958). 


Yuille, I. M., “On the Constraint at the Ends of Ship's 
Seructural Members , Trans. 1. NegA., vol, 94$ pp. 12-40 (1952), 


Vedeler, Georg, "Calculations of Beams", Trans, I, N. A., 
VoOlweo2, Dp. 20-58 (1950). 


Clarkson, J., "On the Transverse Strength of Ships", European 
Snap ptisaine, No. 2, 1958, 

Shanley, F. R., Weight-Strenath Analysis of Aircraft Structures, 
Mee  Yorkssicecraydiitll boon. Co,, Inc,, 1952. 





A Guide for the Analysis of Shin Structures, Edited by Thien Neh, 
U. S,. Department of Commerce, Washington 25, D. C., 1960. 


ash) 























thes 


Ki \) | | 
"| tif | mm 
i ii} Wi | } 
J NA MA i iH i 
Mi 


HN ii mn viauit 
\ HAN i 


HHH | | 
NV 
i} 


il \ 

Wh lV | 1I| 
Wad Ht it 
rs | 

i) 


} 
|| 
i Hn 


\ 
Ni i LAN i Hi 
WALT Witt 


J 4 





